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The  purpose  of  this  study  was  to  desig^n  a supercriti- 
cal airfoil  test  model  which  would  bo  used  to  irivesti{iate 
the  effects  of  cooling  on  boundary  layer  stability  in  sub- 
sonic -flow.  Previous  investigations  have  indicated  that 
transition  to  turbulent  flow  v;as  delayed  by  cooling. 

This  report  is  limited  in  scope  to  the  model  design 
and  test  parameter  specification.  I liopo  that  tliis  v/ork  will 
be  found  complete  and  self-sufficient  by  the  student  v;ho  un- 
dertakes the  actual  testing  of  the  model.  Anyone  who  is  in- 
terested in  determining  the  heat  transfer  for  flew  over  an 
arbitrary  body  with  constant  surface  temperature  v/ill  find 
the  computer  program  in  Appendix  A to  be  useful. 
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Abstract 


A wind  tunnel  test  model  of  a supercritical  airfoil 
v/as  desigiied  to  investigate  the  v/all  cooling  effect  on  sub- 
sonic boundary  layer  stability.  A DGl'iA  523  airfoil  section 
was  employed.  The  model  v/as  designed  to  have  surface  tem])er- 
ature  instrumentation  and  a liquid  nitrogen  cooling  system. 

Heat  transfer,  aerodynamic  loads  and  stresses,  and 
instrumentation  v/ere  analyzed  for  the  proposed  test  condi- 
tions. A computer  program  was  developed  to  analyze  the  forced, 
convective  heat  transfer  over  a two-dimensional  body  v/ith  a 
constant  wall  temperature.  The  program  utilized  an  integral 
method  to  compute  local  Stanton  numbers.  Local  heat  flux  and 
total  heat  flow  were  predicted  for  a f.5ach  number  of  0.7.  Rey- 
nolds numbers  of  0.923  x 10^  and  1.673  x 10^,  and  cooling 
ratios  fi-ora  1.000  to  0.82^1-.  The  stress  analysis  consisted  of 
applying  beam  bending  theory,  along  with  some  simplifying 
assumptions,  to  the  model.  Construction  drawings  and  speci- 
fied test  conditions  for  Mach  numbers  of  0.3.  0.5,  and  0.7 
are  included.  The  proposed  tests  are  to  be  conducted  in  the 
subsonic  test  section  of  the  Trisonic  Test  Facility  at  Wright- 
Patterson  AFB , Ohio. 
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Analysis  and  Design  of  a Cooled  Supercritical 
Airfoil  Test  Model 

I.  Introduction 


Background 

In  1974  Boehraan  and  Mariscalco  of  the  Dayton  Research 
Institute,  University  of  Dayton,  Dayton,  Ohio  began  to  study 
the  problem  of  determining  the  conditions  under  which  flov/ 

in  a cooled  compressible  la.minar  boundary  layer  becomes  un-  ) 

stable  (Ref  2) . This  research  was  under  contract  from  the 
Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson  AFB , 

Ohio.  Using  a computer  program  developed  from  parallel  linear 
stability  theory,  Boehman  looked  at  boundary  layers  in  shock 
tube  induced  flow  and  subsonic  wind  tunnel  flov/. 

This  theoretical  problem  is  still  not  totally  solved. 

A full  description  of  the  stability  problem  requires  the 
solution  of  the  Navier-Stokes  equations  v/ith  boundary  con- 
ditions for  transition.  This  set  of  equations  has  not  been 
solved  to  date.  Even  linearisation  of  these  equations  using 
small  perturbation  theory  yields  equations  that  are  not 
solvable.  Finally,  Boehman  further  assumed  t)iat  the  flow  is 
locally  parallel  and  reduced  tlie  problem  to  a single  ordinary 
differential  equation.  He  found  that  moderate  cooling  re- 
sulted in  a significant  increase  in  stability  of  subsonic 
boundary  layers.  As  a result,  Boeiiman  recommended  that  an 
experimental  test  program  be  developed  to  study  the  effects 
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of  surface  cooling  on  boundary  layer  transition. 

Experimentation  is  needed  to  correlate  with  the  linear 
stability  theory  and  validate  the  assumptions  made  during 
the  development  of  Boehman's  theory.  Very  little  experimental 
work  has  been  performed  on  boundary  layer  stability  in  sub- 
sonic flows.  However,  there  has  been  extensive  work  in  super- 
sonic flow  and  shock  tube  flow  (Ref  l6  and  3)  which  gives 
insight  to  the  stability  problem.  Although  similar  results 
would  be  expected  for  subsonic  flow,  experimentation  is  re- 
quired to  confirm  this  fact. 

The  fact  that  skin  friction  and  the  rate  of  heat  trans- 
fer are  an  order  of  magnitude  greater  in  turbulent  flow  than 
in  laminar  flow  at  the  same  Reynolds  number  makes  the  ability 
to  predict  and  control  the  transition  to  turbulent  flow  of 
great  use  in  the  design  of  aerospace  vehicles.  Most  impor- 
tantly, the  control  of  transition  will  result  in  reduced 
aerodynamic  drag  with  the  following  possible  benefits: 

1 . Increased  aircraft  performance 

2.  Extended  range 

3.  Increased  payload 

4.  Improved  fuel  economy 

5.  Use  of  hydrogen  fuel  feasible  (fuel  cools  wing) 
Problem 

The  problem,  which  was  proposed  by  the  Air  Force  Flight 
Dynamics  Laboratory,  was  to  design,  construct,  and  test  a 
supercritical  airfoil  model  in  the  subsonic  test  section  of 
the  Trisonic  Test  Facility  at  Wright-Patterson  AFB,  Ohio. 
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The  test  v/as  to  invosti{^ate  the  wall  cooling  effect  on  boun- 
dary layer  stability  and  trcinsition . The  purpose  of  the  test 


v/as  to  increase  the  understanding -of  boundary  layer  stabili- 
ty, to  obtain  clues  that  might  advance  the  stabi] ity  theory, 
and  to  correlate  with  the  present  theory. 

.The  model  was  to  have  a DSMA  523  airfoil  section  with 
the  capability  of  being  cooled.  This  White omb-type , super- 
critical airfoil  profile  v/as  developed  by  tlie  McDonnell 
Douglas  Corporation.  The  non-dimensional  surface  coordinates 
of  this  airfoil  were  obtained  from  Ref  8:738.  The  angle  of 
attack  v/as  eliminated  as  a variable  by  specifying  a constant 
zero  angle  for  the  airfoil  section  chord  line.  This  model 
v/as  required  to  have  a constant  airfoil  section,  to  spa.n 
the  entire  24  in.  wind  tunnel,  and  to  be  mounted  on  two  "L" 
shaped  struts.  Major  components  of  the  model  and  struts  are 
illustrated  in  Fig  1 . 

Scope 

The  problem  was  divided  into  tv/o  phases.  The  present 
investigation  completed  the  analysis  and  design  phase  whicii 
included  specification  of  testing  procedures.  The  second 
phase  v/hicli  consists  of  model  construction  and  testing  will 
be  performed  as  a separate  independent  study. 

An  analytical  study  was  required  to  design  the  model 
and  its  cooling  system.  This  analysis  required  the  solution 
for  the  forced,  convective  heat  transfer  in  subsonic  flow 
over  the  two-dimensional  wing.  A computer  progi-am  was  devel- 
oped from  the  reference  temperaLure  approacl)  to  variable 
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Fig.  1 Concept  of  the  Model  v/ith  Moiinting  Struts 

property  flow  which  used  integral  equations  to  determine 
the  local,  laminar  and  turbulent  Stanton  numbers  for  constant 
wall  temperature  flow  over  aji  arbitrai-y  body  (Ref  9:226  and 
247).  Estimates  of  heab  conduction  in  the  model  were  also 
necessax-y.  Tlie  )ieat  transfer  axialys^s  which  v/a.^  performed 
for  only  the  most  stringent  tost  conditions  appears  in  detail 
in  Section  II. 

The  model  design  required  determination  of  t]io  forces 
and  moments  on  the  iriodel,  a stress  analysis,  and  construction 
drawings.  The  stress  analysis  made  use  of  simple  beam  bending 
theory.  Details  of  the  stress  analysis  which  was  performed 
only  for  the  most  stringent  test  conditions  ajqioar  in  Section 
III. 
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The  teat  was  designed  to  measure  the  movement  of  boun- 
dary layer  transition  as  a function  of  wall  cooling  i-elative 
to  the  no  cooling  condition.  The  stability  results  of  Boeh- 


i 
I 

[ man  (Ref  2:5?)  were  used  in  the  forced,  convective  heat 

I transfer  analysis.  The  test  conditions  which  are  specified 

I in  Section  IV  were  limited  to  combinations  of  the  following 

\ variables  for  the  ranges  indicated; 

1.  Reynolds  number  (1  .673  x 10^  and  0.923  x 

1 2.  Mach  number  (0. 7-0.3) 

! 

i 3-  Cooling  ratios  (1.000  - 0.824) 
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I I . Heat  Transfer  Analysis 

Boimdar y Layer  Force d Convection 

General  Assumptions.  Simplifying  assumptions  were 
made  to  the  model  so  that  existing  solutions  to  the  heat 
transfer  ])roblem  could  he  used.  With  the  proposed  model 
description,  the  fluid  flow  v/as  assumed  to  be  tv/o-dimensional . 
It  was  assumed  that  only  the  most  stringent  test  condition 
of  M=0.7  needed  to  be  analyzed  since  aerodynamic  heating 
and  heat  flux  increase  as  a function  of  Mach  number.  A con- 
stant wall  temperature  v/as  specified  since  it  eliminated 
tempera tu.re  gradients  as  an  influencing  factor  iji  boundary 
layer  transition. 

Boehman  reported  i-'esults  that  give  the  transition 
Reynolds  number  as  a function  of  the  wall  cooling  ratio, 

, for  a family  of  subsonic  Mach  numbers  (Ref  Zt  57)  • 

The  results  apply  to  a flat  plate.  It  was  assumed  that  these 
results  (Fig  1)  could  be  used  as  a first  approximatioii  to 
the  transition  on  the  model. 

At  tlie  .stagnation  point  on  the  airfoil,  the  equations 
used  were  singular.  In  order  to  compute  the  convective  con- 
ductance, h,  the  flow  was  assumed  to  be  approximated  by  the 
potential  flov/  over  a cylinder  (Eq  1). 

Ug  = 2.U^x/R  Cx/R<<1-0) 

For  the  range  of  Mach  numbers  of  interest,  M“0.3  to  0.7, 
compressibility  could  not  be  neglected.  Thus,  the  flow  was 
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CRITICAL  REYNOLDS  NO.,  ('^crit'^'^’criT  ^ 


Fig.  2 Relative  Stability  of  Boundary  Layers  as  a Function 
of  Critical  Reynolds  Number  and  Cooling  Ratio 
(From  Ref  2:57) 

considered  to  have  variable  properties.  A reference  tempera- 
ture approach  was  used  in  which  all  projierties  v/ere  computed 
locally  using  Eckert's  reference  temperature, 

/ = + 0.  Z ?.  ( 2 ) 

Theory . The  solution  to  the  convective  heat  ti'ansfer 
analysis  could  have  been  obtained  by  solving  the  bouiidary 
layer  momentum  equation  which  establishes  the  velocity  field 
and  tlien  solving  the  boundary  layer  energy  equation.  This 
method  v/as  complicated,  and  tlie  accuracy  of  results  required 
did  not  merit  its  use.  Simpler  metliods  were  available  when 
the  proper  approximations  were  made  to  the  problem.  Use  of 
the  simpler  methods  required  that  tlie  flow  over  the  model  be 
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analyzed  in  three  ree,iona,  utagnation,  laminar,  and  turbulent 
flow . 


Kaye  (Ref  9:223)  dovclo])ed  a method  of  computing  the 
local  Stanton  number  for  laminar  flow  over  a constant  tem- 
perature body  of  arbitrary  shape . Tliis  method  v/as  based  on 
the  assumption  that  the  variation  of  conduction  thickness, 
for  flow  over  an  ai’bitrary  shape  \.'as  of  the  same  form 
as  the  similarity  solutions  for  Vv^edge  flow.  The  function, 
which  expi-essed  the  variation  of  conduction  thickness  in 
air,  Pr=  O.7,  for  v/odge  flow,  was  defined  numerically.  Smith 
and  Spalding  (Ref  13:60)  approximated  this  function  by  a 
straight  line  v/hich  fits  the  v/edge  solutions  exactly  for 
the  stagnation  point  and  the  flat  plate.  The  resulting  ex- 
pression wtiS  solved  for  the  conductioj)  thickness.  Using  the 
definition  of  conduction  thiclniess  and  Stanton  number,  the 
following  expression  was  derived  for  the  local  Stanton  num- 
ber; 


. , „ 0-5  / . \ O.'rSS 

0.4 1 & / 


O.K 


(3) 


I For  the  turbulent  case  of  flow  over  a constant  temper- 

t- 

ature  body  of  arbitrai'y  shape,  Ambrok  (Ref  1:1979)  suggested 
I a solution  of  the  integral  energy  equation  without  h£iving 

to  solve  the  momentum  equation.  This  method  was  developed 
I by  Kays  (Kef  9:245-2^7).  Kays  started  with  the  solution  to 

I heat  convection  for  the  constant  temper-ature  flat  plate  in 

I 

turbulent  flow.  This  solution  whicli  was  expressed  in  terms 


I 


of  tho  Stanton  nurabor  was  used  in  the  energy  integral  equa- 
tion for  a flat  plate  to  develop  an  expression  for  the  en- 
thalpy thickness,  A2  • The  resulting  relationsliip  betv/oen 
X and  A 2 used  to  develop  a Stanton  number  v/hich  v/as  a 
function  of  local  parameters  only.  Tliis  Stanton  number  re- 
lationship was  assumed  to  hold  for  any  turbulent  flow.  The 

general  energy  equation  v.-as  integrated  with  tlie  heat  flux, 

>11  . 

q^,  in  the  equation  expressed  in  terms  of  the  convective 
conductance  which  was  found  from  the  assumed  Stanton  number 
relationship.  The  resulting  expression  v;as  soli'-ed  for  the 
following  local  Stanton  number: 


St.  =0.0^95 


[r 


e 


u. 


dx  ] 


o.a 


(^) 


At  the  stagnation  point,  Eq  3 for  the  local  Stanton 
number  did  not  apply  since  Ug  = 0.  However,  Eq  3 v/as  used  to 
solve  for  the  convective  conductance,  h,  at  x=0.  After  sub- 
stituting the  assumed  velocity  relationship  from  Eq  1 , the 
expression  for  h was  integrated  witVi  the  following  results: 


,001 


c,. 


R 


The  value  of  h was  now  independent  of  the  values,  x=0  and 

u^=0 . 
e 


In  ali  three  regions  of  flow,  the  local  heat  flux  v/as 
determined  from  the  convective  conductance  obtained  from  Eq 
3.  and  5 !ind  the  knov/n  local  temperature  difference  by 


(6) 


C,  = h,  ( T,  - ) 

Comnutcr  Prop^ran  Iioqic.  A computer  program  call  CON- 
HEAT  was  written  to  analyze  the  forced,  convective  heat 
transfer  0]i  tlie  airfoil  model  . This  program  utilized  Eq  2 
th’^ough  Eq  6.  The  program  was  designed  to  he  general  in 
nature  and  could  be  utilized  to  compute  the  convective  heat 
trajisfer  for  any  tv/o-dimensional  body  with  a constant  surface 
temperature . Any  further  restrictions  on  the  application  of 
the  program  will  be  noted  during  the  discussion  of  program 
logic.  A detailed  description  of  the  input/ output  format  is 
given  in  Appendix  A along  with  a complete  listing  of  the 
program . 

Information  required  to  solve  the  convection  problem 
included  freestream  conditions,  a pressure  distribution  over 
the  model,  model  geometry,  and  the  desired  cooling  ratio. 

The  required  freestream  conditions  consisted  of  the  velocity, 
prcssui'e,  temperature,  and  viscosity.  Pressure  coefficients 
were  required  on  the  surface  at  intervals  of  5 percent  of 
the  chord  starting  at  the  stagnation  point.  This  resulted  in 
21  points  on  the  upper  surface  and  21  points  on  the  lower 
surface  being  defined  for  use  in  determining  tlie  total  rate 
of  heat  transfer.  A surface  coordinate  system  was  defined, 
and  the  distance  between  pressure  coefficients  on  the  sur- 
face was  determined.  Additional  geometry  required  included 
the  chord  length,  span,  and  leading  edge  radius.  The  cooling 


I 
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ratio,  , rjpecified  the  surface  tenixicrature  i'or  the  giv- 

en freestreain  temperature,  and  the  critical  Reynolds  number 
was  obtained  from  Fig  2. 

The  first  section  of  the  program  computed  the  proper- 
ties at  the  discrete  locations  of  the  pressure  coefficients 
for  the  upper  cind  lovvei-  sui-faces  io  addition  to  tlie  free- 
stream  density  and  the  recovery  factors.  The  laminar  recovery 
factor  was  defined  by 

r = ^/P^  ( 7 ) 

Lam 

The  turbulent  recovery  factor  was  defined  by 

r 

Tar 

The  freestream  density  was  computed  from  th'e  perfect  gas  lav/. 

The  velocity,  temperature,  and  pressure  at  each  point  were 
determined  by  using  inviscid  relations  wiiich  did  not  include 
the  effect  of  boundary  layer  displacement  thickness . 

The  second  section  computed  the  Eckert  reference  temp- 
erature and  the  properties  based  on  this  temperature  for  the 
upper  surface . The  adiabatic  wall  temperature  which  is  used 
to  compute  Eckert's  temperature  was  found  by 

T.w=  t 

whei’o  r was  the  value  fox'  laminar  flow  from  Eq  7*  Then  Eck- 
ert's temperature  was  computed  from  Eq  2.  The  viscosity  was  j 

computed  using  Eckert's  temperature  in  Sutherland's  equation. 
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T^.  -t-is?,./;') 
T*+i3^.V;:/ 


The  local  donalty  v/as  computed  using  the  pox-fect  gas  law. 

A local  Reynolds  number  was  computed  and  compax-ed  with  tlio 
critical  Reynolds  number.  If  tlie  value  of  tlie  critical  Rey- 
nolds number  had  been  exceeded,  then  the  computations  v;ero 
repeated  with  the  value  of  r defined  by  Eq  8 for  the  remain- 
ing points  on  the  surface  in  turbulent  flow.  Fox’  the  third 
section,  the  same  pi'-ocedures  as  in  the  second  section  were 
repeated  fox-  the  lov/er  surface. 

The  fourth  section  computed  the  local  Stanton  number, 
local  convective  conductance,  and  local  heat  flxix  at  each 
point  on  the  upper  surface.  Equation  5 was  used  to  compute 
the  convective  conductance  at  the  stagnation  point.  Equation 
3 was  numei'-ically  integrated  by  the  trapesoidal  x''ule  to  ob- 
tain the  laminar-  Stanton  number.  After-  the  critical  Reynolds 
number  had  been  exceeded,  Eq  4 v/as  used  instead  of  Eq  3 to 
compute  the  local  turbulent  Staixton  number.  The  follov/ing 
relation  was  used  to  compute  the  convective  conductance  fi'-om 
the  Stanton  number: 

= St ^ ^ Cp  (11 ) 

The  local  heat  flux  was  determined  from  Eq  6.  Fox'-  tiie  fifth 
section,  the  same  procedures  as  in  the  fourth  section  wore 
repeated  for  the  lower  sui-'face . 

The  last  secti.on  integrated  the  local  heat  flux  by  the 
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trapeisoidal  rule  for  the  upper  and  lower  surfaces.  These 
results  were  added  for  both  surfaces  with  the  final  result 
being  the  total  rate  of  heat  transfer.  Other  results  which 
were  part  of  the  program  output  v/ere  the  local  conditions  at 
each  point  and  a plot  of  the  heat  flux  distribution.  Figure 
3 summarizes  the  program  logic. 

Design  Data.  The  most  stringent  test  condition  deter- 
mined the  design  data  which  was  used.  A Mach  num.ber  of  0.7 
was  selected.  This  corresponded  to  a velocity  of  8l3  ft/sec 
at  the  wind  tunnel  temperature  of  100  F.  The  freostream  vj.s- 
cosity  from  Eq  10  v;as  0.37049  x 1 0~°lb^-sec/ft  . The  free- 
stream  pressure  was  determined  from  a method  described  in 
Section  IV.  This  pressure  and  the  corresponding  Reynolds 
number  wero  1320  lb^,/ft^  and  1 .673  x 10^  respectively.  The 
desired  cooling  ratio  was  0.620.  The  airfoil  geometry  v.-as 
obtained  graphically  from  a full  scale  drawing  of  the  airfoil 
section.  The  last  item  required  v/as  the  pressure  distribution 
which  was  obtained  from  Ref  8:739* 

Results . A plot  of  the  heat  flux  for  Re=1.673  x 10^ 
and  T^T^  =0.620  appears  in  Fig  4.  The  negative  values  of 
the  heat  flux  indicated  heat  transfer  from  the  boundary  layer 
to  the  airfoil.  There  was  a large  flux  near  the  leading  edge 
which  v;as  expected  for  the  stagnation  region.  The  flux  dropped 
off  quickly  in  the  laminar  flow  region.  The  jump  in  heat  flux 
at  x/c=0.25  indicated  transition  to  turbulent  flow.  The  tur- 
bulent heat  flux  which  was  larger  than  the  laminar'  flux  grad- 
ually decreased  to  the  trailing  edge  of  the  airfoil. 
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INPUT  (Appendix  A): 

FREESTREAM  GOND IT iONE 
COOLING  RATIO  DESIRED 
AIRFOIL  GEOI.IETRY 
PRESSURE  DISTRIBUTION 


i 

COI.IPUTE  SURFACE  CONDITIONS 

1 

1 

1 

(Eq  2,  7 - 10). 

* 

'^e'  ^e’  ^ ^ 

USING  ECKERT'S  TErJP , 

T’ 

COPIPUTE  LOCAL  STANTON  NO.: 

LAMINAR  (Eq  3) 

u - - -MLS  A (p  M _ 

o‘u  c“  r c"/  * I 

TURBULENT  (Eq  4) 

— 0 4 ■**  0. 2_ 

<-+  = -O.OZOE  ?r  _ M _ 

''  Li*  ? 


Fig.  3 Flov/  Dififiraiii  for  GONiiEAT  Pi-o^qrarn 


INTEGRATE  KEAT  FLUX  FOR  TOTAL  liEAT  FLOWi 


QrC  .»  ,1 

= ) c\«(x)ax 


OUTPUT  (Appendix  A): 


Re  (LOCAL) 


Fii".  3 (continued)  Flow  Diafjram  for  CONIUiAT 


1.5 


Further  results  of  the  convection  analysis  appear  in 
Table  I.  This  table  describes  the  local  conditions  for  the 
uppei'  and  lower  surfaces  as  a function  of  location  in  the 
surface  coordinate  system.  The  convective  conductance  and 
Stanton  number  followed  the  same  trend  as  was  noted  for  the 
heat  flux.  As  a result  of  integrating  the  local  heat  flux 
over  tlie  airfoil  surface,  the  total  heat  transfer  rate  v/as 
17-03  B/sec . 

Heat  Conduction  Estimates 

Heat  conduction  estimates  v/ere  performed  to  aid  in  the 
proper  selection  of  construction  materials  for  the  airfoil 
surface  and  insulating  end  sections  and  to  determine  the  pro- 
per thicknesses  for  these  materials.  A detailed  analysis  of 
the  steady  state  heat  conduction  problem  would  have  required 
the  solution  of  the  two-dimensional  Laplace's  equation, 

for  a constant  v/all  temperature  and  specified  heat  flux  as 
boundary  conditions.  The  distribution  of  heat  flux  varies  as 
a function  of  freestream  conditions  ar.d  cooling  ratio  as 
can  be  seen  in  Appendix  B.  Tims,  after  a solution  v/as  ob- 
tained, it  v/ould  be  valid  for  only  one  test  condition.  Using 
this  solution  to  design  the  constant  temperature  airfoil  sur- 
' face  would  result  in  a model  with  a constant  temperature 

surface  only  at  that  test  condition.  In  addition,  model  con- 
struction would  be  complicated  by  the  contoured  surfaces 
required  to  achieve  the  constant  tempcn’ature . 
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R^RFOI^  HEf^7  nyx  (ajU/SEp-SQF} 


X/C  (LOCflTION/CHORD) 


Fig.  4 Airfoil  Heat  Flux  at  the  Doai^-ji  Condition 
1V1=0.7,  Ro=1.673  X 10^,  ^^0.620 


17 


Table  I 


Airfoil  Local  Properties  at  the  Design  Condition 
iVl=0.7,  He=1.673  X 10^,  T yT^  =0.620 


As  an  obvious  concussion  to  the  method  discussed,  a 
simple  estimate  v/as  applied.  The  Fourier  equation  for  heat 
conduction  v/as  approximated  by 

• “ Ic 

(12) 

where.!  is  the  thickness  normal  to  the  airfoil  surface. 
Equation  12  v/as  used  with  the  results  of  the  convective  heat 
transfer  analysis.  An  average  heat  flux  v/as  found  for  the 
upper  airfoil  surface  only  since  this  v/as  the  surface  of 
primary  interest.  This  heat  flux  v/as  computed  by  taking  the 
mathematical  average  of  the  flux  at  each  point  excluding  the 
stagnation  point  from  Table  I. 

For  the  case  of  the  airfoil,  the  maximum  acceptable 
temperature  difference  through  the  surface  was  specified  to 
be  3 F . The  thickness,  t,  v/as  assumed.  Then  a temperature 
difference  v/as  computed  using  the  average  heat  flux.  This 
tempei'ature  difference  was  compared  to  the  limit  to  determine 
acceptability . 

The  minimum  insulating  end  section  thickness  for  the 
airfoil  which  was  required  to  prevent  the  w.ind  tunnel  v/indov/ 
from  breaking  or  condensat.ing  moisture  was  determined  by  a 
heat  flow  balance.  The  cooling  of  the  windov/  caused  by  con- 
duction in  the  model  was  equated  to  the  heat  convection  to 
the  v/ind  tunnel  window.  The  flow  over  the  window  was  assumed 
to  be  fully  developed  turbulent  duct  flow.  A temperature  of 
70  F was  required  to  be  maintained  in  the  window  glass.  T])is 
heat  balance  was  solved  for  the  required  thickness. 
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I Selections  of  construction  materials  and  material  thick- 

; nesses  v/ere  made  using  these  estimates.  The  average  heat 

\ p 

i flux  on  the  upper  surface  of  the  airfoil  v/as  5 •^■*■51  B/sec-ft  . 

t 

f Aluminum,  which  has  a thermal  conductivity  of  O.O3806  B/sec- 

r 2 

! ft  -R/ft,  was  selected  for  the  airfoil  center  section.  A 

L 

I 0.187 -in.  thick  surface  of  aluminum  caused  a temperature 

I difference  of  2.23  F when  computed  from  Eq  12.  Thus,  the 

j airfoil  was  designed  with  a constant  thickness  of  0.18?  in. 

White  oak  v/hich  has  a thermal  conductivity  of  1 .38  x 10  ^ 

I 2 

I B/sec-ft  -R/ft  was  selected  for  the  insulating  end  sections. 

A 0.0122  in.  minimum  thickness  is  required  to  keep  the  v/in- 
I dow  at  70  F.  However,  due  to  structural  considerations , the 

end  sections  v/ere  designed  to  he  0.81 3 in.  thick. 

1 

: ! Cooling  System  Design 

The  airfoil  cooling  system  v/as  designed  to  permit  var- 

I iation  of  the  cooling  ratio  during  the  wind  tunnel  test.  The 

! 

j system  had  to  be  capable  of  cooling  the  model  to  -II3  F which 

I 

i corresponds  to  the  cooling  ratio  of  0.620  and  had  to  minimize 

i temperature  gradients  in  both  the  spanv/ise  and  chordv/ise  di- 

rections . 

I Liquid  nitrogen  was  selected  as  the  coolant.  This  se- 

i 

j lection  was  based  on  the  following  properties  (Ref  4:7^)  at 

i ' one  atmosphere  of  pressure: 

Tb,|=-  320  F 

i 

' ^ = 50.-)  6 
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The  cooling  system  was  designed  to  operate  on  liquid 
nitrogen  only.  The  mass  flow  rate  of  liquid  required  to  cool 
the  airfoil  at  the  design  condition  was  estimated  by 


m=  (13) 

where -Q  is  the  total  heat  transfer  rate  from  Table  I. 

The  inlet  area  required  to  carry  the  liquid  v/as  deter- 
mined by  assuming  that  the  density  was  constant  and  that  the 
average  velocity  of  the  fluid  v/as  20  ft/sec.  To  ensure  a 
safety  factor  in  the  cooling  capacity,  this  inlet  area  v/as 
increased  by  120  percent. 

The  airfoil  is  cooled  by  spraying  liquid  nitrogen  on 
the  inner  walls  of  the  model  from  nozzles.  Most  of  the  energy 
is  absorbed  by  vaporization,  and  the  resulting  gas  is  expelled 
into  the  wind  tunnel.  The  most  complicated  part  of  the  cool- 
ing system  design  was  locating  and  sizing  these  nozzles. 

The  cooling  system  design  was  based  on  the  predicted 
total  heat  transfer  rate  and  the  heat  flux  distribution  at 
the  design  condition.  The  total  heat  flew  of  17.03  B/scc  was 
used  to  approximate  the  mass  flow  of  liquid  nitrogen  required. 
From  Eq  13»  this  mass  flow  rate  was  O.I98  Ib^^j/sec.  The  inlet 
area  to  the  airfoil  which  v/as  computed  and  multiplied  by  the 
safety  factor  was  0.0622  in  . The  nozzles  which  are  pipes  of 
1/8  in.  nominal  size  were  located  as  shown  in  Fig  5*  This 
location  was  based  on  the  heat  flux  distribution  in  Fig  4. 

The  required  inlet  area  was  achieved  by  drilling  the  nozzles 
witli  a number  76  bit  for  a total  of  198  holes.  Tiie  holes  were 
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Lov;er  Section  with  Cooling  Nozzles 
foil  Boss 


distributed  over  tlie  entire  length  of  the  nozzles  v/ith  a 
larger  number  of  holes  in  areas  of  high  heat  flux.  Forty- 
four  percent  of  the  holes  must  be  in  the  leading  edge  nozzle. 
The  middle  nozzle  must  contain  3^  percent  of  the  holes.  Tv/en- 

ty-two  percent  of  the  holes  must  be  in  the  trailing  edge 

i 2 

nozzle.  The  required  exhaust  area  v/as  0.864B  in  , which  is 

! 

j times  the  area  of  the  inlet  to  ensure  unrestricted  flow. 

I The  remainder  of  the  system  consisted  of  the  liquid 

I 

: nitrogen  supply  source  and  connections  through  the  mounting 

struts  of  the  model.  One  strut  contained  the  liquid  nitrogen 

supply  line.  Both  struts  contained  exhaust  pipes  for  the 

nitrogen  gas  as  illustrated  in  Fig  1 . The  nitrogen  supply 

for  the  system  consisted  of  a dewar  of  liquid  nitrogen  with 

insulation,  a compressed  nitrogen  gas  bottle  at  2250  psi, 

and  a pressure  regulator.  Regulated  compressed  gas  pressur- 
i 

I ized  the  dewar  to  permit  adjustment  of  the  mass  flov/  r3.te . 

i Pressure  relief  values  were  used  to  prevent  ovei'  pressuriz- 

I ing  the  model,  piping,  and  dewar. 
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III.  A.ii'foil  Test  Motlel  DeGi/’:a 


Aerodynamic  Loadfj 

Test  loads  on  the  airfoil  and  the  mounting  struts  v/ere 
predicted.  Lift,  drag,  and  pitching  moment  act  on  the  airfoil. 
The  only  force  acting  on  the  syinmetrical  struts  is  drag.  A 
small 'number  of  known  experimental  results  for  the  DSIM  523 
airfoil  served  as  initial  conditions  for  the  computation. 

These  known  values,  which  were  obtained  from  Ref  8,  are  the 
following: 

M = 0.82 

Re^  = 3 X 10^ 

= 0.51 
= 0.015 

G , Pressure  distribution 
P 

The  angle  of  attack  and  moment  coefficient  for  the  conditions 
were  unknown.  The  problem  was  to  find  the  force  and  moment 
coefficients  for  M=0.?  and  Re  =1.673  x 10^  at  the  zero  angle 
of  attack . 

Angle  of  Attack.  Equations  were  derived  which  relate 
the  force  and  moment  coefficients  to  the  angle  of  attack  and 
were  used  to  determine  the  angle  of  attack  for  the  known 
experimental  conditions.  The  forces  and  moments  that  act  on 
an  airfoil  are  depicted  in  Fig  6. 

The  summation  of  forces  in  the  x-  and  z-  directions  and 
the  summation  of  tlie  moments  about  the  leading  edge  must  be 
equal  to  zero  to  satisfy  equilibrium.  A summation  of  forces 
in  the  x-  and  z--dii.'cction  using  Fig  6 yields; 


Forces  and  Moments  on  a Typical  Aii'foil  Section 


c.osc<. 


about  the  leading  edge  gives 


Finally,  summation  of  moment 


Equation  wa: 


olved  for  the  angle  of  atta  k to 


obtain 


The  angle  of  attack  was  determined  to  be  '->< 


wcT-e  known. 

Detormlnation  of  Quarter  Chord  Fitchin/^  Moment  Coeffi - 
cient . Using  the  plot  of  jjrossuro  distribution  of  C^  = 0.51i 
the  moment  coefficient  about  the  quarter  chord  point  was  de- 
termined to  be  ^,y^=0.123.  Details  of  the  method  used  can 
be  found  in  Ref  7:83-85  and  consisted  of  numerically  inte- 
grating the  pressure  distribution  for  directions  both  normal 
and  perpendicular  to  the  chord.  Then,  an  appropriate  coordi- 
nate transformation  v/as  performed.  This  result  was  verified 
by  the  solution  of  Eq  l6  w’ith  very  close  agreement. 

Determination  of  Theoretical  Coefficients.  The  prev- 
iously derived  equations  and  results  were  used  to  find  the 
theoretical  coefficients  at  the  zero  angle  of  attack.  It  was 
assumed  that  the  drag  coefficient  which  v/as  presented  contain- 
ed only  pressure  drag.  Thus,  it  was  necessary  to  a.dd  s]:in 
friction  drag.  The  skin  friction  coefficient  was  calculated 
using  the  method  of  Van  Driest  (Ref  15)  which  computes  t!ie 
comipressible , turbulent  skin  friction  coefficient  for  a flat 
plate  as  a function  of  heat  transfer.  The  adiabatic  condition 
was  assumed.  The  resulting  total  drag  coefficient  was  0.03^'-. 

The  ideal  lift  curve  slope  for  a flat  plate  v/as  assumed 
to  apply,  and  the  following  relation  was  derived: 

v/hero  "i"  denotes  initial  conditions.  Equations  l4,  16,  and 
l8  define  tlie  relationship  between  the  coefficients  and 
angle  of  attack . 
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Determination  of  Fox-cen  and  Momont . T)ie  forces  and 
moment  on  tlie  airfoil  could  be  comy^uted  v/hen  the  dynamic 
pressure  and  the  wing  plant orm  area  v/ere  specified.  The  dy- 
namic pressure  was  expressed  as 


where  =1  .A  for  air  and  M,  and  P v/ere  defined  for  the 

design  condition.  The  airfoil  chord  was  specified  to  be  0.750 

ft  and  the  span  was  specified  to  be  2.0  ft.  Thus,  the  plan- 

2 

form  area  was  1 .50  ft  . 

Strut  Drag.  Only  tho  vertical  portion  of  the  strut  was 
considered  since  the  horizontal  portion  in  tho  v/ake  of  the 
vertical  strut  was  assumed  to  contribute  negligibly  to  the 
total  drag.  Since  the  sxruts  were  symmetric  and  pai^allel  to 
the  flow,  the  only  force  that  occurred  v/as  drag.  The  Reynolds 
number  v/as  487,800  for  the  2.625  in.  choi-d  of  the  strut.  The 
thickness  ratio,  t/c , was  0.57  or  57  percent.  Based  on  a 
critical  Reynolds  number  of  110,000  for  an  approximately  50 
percent  tlriok  bluff  body  as  suggested  by  Hoerner  (Ref  6:6-5), 
t}29  location  of  transition  was  0.0493  fi--  The  laminar  and 
turbulent  skin  friction  coefficients  wei-e  computed  by 


=a074/R.’f 


A weighted  mean  value  based  on  length  v/as  found  for  the  skin 
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friction  coefficient.  The  totci.l  drag  coefficient  v/as  esti- 
mated by  (Ref  6:6-9) 

Cj  jC^  = 4+2  % + l20(+/p,j'"  (20) 

where  is  based  on  frontal  ai-ea  and  t/c  is  the  thickness  to 
chord  -ratio. 

Results . At  M=0.7  and  o<:  = 0°,  the  aerodynamic  loads 
which  were  computed  from  the  theoretical  coefficients  for 
the  airfoil  are  the  following: 

L - 299.0  Ibj:. 

D = 12.9  Ib^ 

M = 56.0  ft-lb^ 

The  total  strut  drag  was  l4.5  Ib^.  The  theoretical  design 
coefficients  for  lift  cind  drag  were  0.440  and  0.019»  respec- 
tively . 

Stress  Ana-lysis 

A stress  analysis  v/as  performed  on  the  entire  airfoil 
model  and  the  wind  tunnel  mounting  struts  to  ensure  that  the 
model  v/as  capable  of  sustaining  the  loads  that  were  predicted 
for  the  design  condition.  Specifically,  the  safety  factor 
was  required  to  be  a minimum  of  4.  Simplifying  assumptions 
reduced  the  analysis  to  an  application  of  beam  bending  theory. 
The  model  and  mount  were  divided  into  major  componorils  (Fig 
1);  the  airfoil,  airfoil  mounting  boss,  v/ind  tunnel  mounting 
struts,  and  strut  mounting  plates.  Th(i  procedures  used  were 
applied  to  all  components. 


28 


Assumptions . Duo  to  the  tv/o-dimensional  flow,  the 
spaiwise  force  and  moment  distributions  on  the  airfoil  model  I 

were  uniform.  The  chordv/ise  force  distribution  was  assumed 
to  be  such  that  the  net  result  acted  at  the  quarter  chord 
point  of  each  section.  This  assumption  is  normal  for  sub- 
sonic ‘flov/.  The  reaction  of  the  mounting  struts  to  the  air 
loads  was  assumed  to  occur  at  the  center  of  the  strut  exhaust 
pipe  and  to  act  as  a point  force  concurrent  with  the  lift 
force.  The  airfoil  v/as  assumed  to  act  as  a simple  beam  with  ' 

a uniform  load  where  the  struts  supported  no  spanv/ise  moments. 

To  find  the  torsional  stress  due  to  the  aerodynamic  moment, 
the  airfoil  was  assum.cd  to  be  an  elliptical  shell  v/ith  the 
major  axis  equal  to  the  chord  and  the  minor  axis  equal  to 
the  maximum  thickness  of  the  airfoil.  The  cross-sectional 
area  and  moment  of  inertia  about  the  axis  normal  to  the  chord 
were  forced  to  be  the  same  as  for  the  airfoil. 

Due  to  model  symmetry,  the  loads  were  divided  equally 
to  each  strut.  Since  the  aii^foil  and  wall  acted  as  end  plates, 
the  flow’  over  the  strut  was  assumec  to  be  two-dimensional  and 
to  create  a uniform  drag  force.  The  forces  and  moment  on  the 
strut  were  assumed  to  occur  at  the  center  of  the  exliaust 
pipe  w'hich  is  located  at  28  percent  chord.  The  analysis  of 
the  strut  v/as  divided  into  the  vertical  portion  and  horizontal 
portion,  where  both  were  assumed  to  act  as  cantilever  beams 
with  concentrated  loads.  This  assumption  was  valid  if  deflec- 
tions at  the  fixed  ends  were  negligible.  All  loads  were  as- 
sumed to  be  carried  by  the  exhaust  pipe  for  tlie  vertical 
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portion  and  by  the  shape  labeled  Section  A-A  in  Fig  7 for 
the  horizontal  portion. 

The  airfoil  mounting  boss  and  boss  plate  were  desigriod 
to  connect  the  airfoil  to  the  strut  and  are  illustrated  in 
Fig  8 and  Fig  9.  respectively.  Any  vertical  transl.ation  of 
the  drag  force  was  assumed  negligible  since  this  force  is 
small  and  the  distance  of  translation  v/as  less  than  one  inch. 
An  equivalent  couple  acting  on  the  centerline  was  assumed  to 
replace  the  aerodynamic  moment. 

The  strut  mounting  plate  was  designed  to  fasten  the 
strut  to  the  wind  tunnel  wall.  The  forces  that  acted  on  this 
plate  were  assumed  to  behave  as  point  forces  acting  at  the 
centroid  of  the  four  mounting  bolts  which  v/as  0.375  in.  from 
the  location  used  to  compute  the  stresses  on  the  horizontal 
portion  of  the  strut.  The  forces  and  moment  on  the  horizontal 
strut  were  assumed  to  act  on  the  plate  since  translation  over 
the  small  distance  v/as  negligible. 

General  Procedures.  Moments  of  inertia  and  centroids 
of  area  w(  re  computed  for  the  cross-section  of  each  component 
from  a method  of  built-up  elements.  Any  complicated  section 
shape  was  divided  into  simple  elemental  areas.  Then,  the 
following  formulas  were  used: 


N 


A = Z 0.; 

csl 

(21) 

N 

(22) 
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Kounting  Strut 
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If  the  moments  of  inertia  were  found  about  a set  of  axes  at 
any  point  other  than  the  centroid  of  the  cross-section,  then 
the  parallel  axis  theorem  v/as  used  to  obtain  the  moments  of 
inertia  about  the  centroidal  axes. 

The  material  selection  was  based  on  the  thermal  and 
structural  properties.  If  after  the  stress  analysis  the  safe- 
ty factor  was  less  than  4,  a nev/  material  with  greater  stren- 
gth was  selected.  The  particular  structural  properties  of 
interest  in  the  analysis  v/ere  the  yield  strengths  and  ulti- 
mate strengths  under  tension,  compression,  and  shear.  These 
were  used  to  compare  with  the  predicted  design  loads  and  to 
define  the  safety  factor. 

Each  component  was  isolated  during  the  analysis.  An 
equivalent  force  system  had  to  be  described  which  produced 
the  same  effect  on  that  component  as  the  reaction  to  the  adja- 
cent component.  The  equivaleiit  system  was  obtained  by  moving 
forces  along  their  lines  of  action  or  by  parallel  translation 
w?iich  required  the  introduction  of  a couple.  Some  uniformly 
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distributed  forces  wore  replaced  by  a point  force  acting  at 
tlie  center  of  the  distributed  load  area. 

Each  component  was  examined  to  define  critical  points 
which  v/ere  suspected  of  sustaining  large  stress  concentra- 
tions due  to  changes  in  cross-sectional  area,  curvature,  or 
distance  from  the  applied  load.  The  critical  points  were 
the  only  locations  where  the  maximum  stresses  and  safety 
factors  v.ere  computed. 

Special  considerations  were  necessary  to  compute  the 
stress  on  the  airfoil  and  the  strut.  The  shear  stress  due  to 
torsion  on  the  airfoil  was  computed  by  a method  detailed  in 
Ref  17; 53*  This  shear  stress  is  given  by 


it  [ cT+''i)’"-Ti 


where 


Q - 0--Q-I  _ -J?l 

- b, 


and  the  subscript,  1 , denotes  the  axial  length  measured  to 
the  inside  of  the  shell.  The  insulating  end  sections  of  the 
airfoil  v/hich  v/ere  assum.ed  to  be  a solid  elliptical  shaft 
required  the  use  of 


[ior  (X  > b) 


The  mounting  struts  had  a sharp  curvature  which  caused  stress 
concentration.  Oberg  (Ref  11:^33)  presented  a method  to  cor- 
rect for  curvature  effects.  This  correction  factor  was  a 
function  of  the  radius  of  curvature. 
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The  compojjGut  v/as  modeled  as  a simple  beam  of  constant 
cross-section.  Vidosic  (Kef  17:31-3^1)  had  compiled  a table  of 
solutions  for  the  i^eactions,  moments,  and  shear  at  any  point 
on  a uniform  beam  for  a variety  of  conditions.  A solution 
v/as  selected  from  this  table  whicii  matched  the  modeled  com- 
ponent- and  equivalent  force  system. 

To  find  the  normal  stress  due  to  the  bending  moment, 
the  flexure  formula, 

v-Mc/I  (29) 

was  used.  The  stress  due  to  each  of  the  forces  or  moment 
was  computed  separately. 

The  resulting  normal  stresses  fi'om  each  force  and  mo- 
ment at  a point  were  combined  by  vector  addition.  Shear 
stresses  were  combined  in  tl-ie  same  manner.  The  remaining 
normal  and  shear  stresses  were  combined  to  find  the  maximum 
stresses  at  the  point.  The  formulas  used  from  Kef  11:4’21  to 
combine  the  stresses  v/ere  the  following: 


(30) 

(31) 

The  maximum  stresses  were  compared  to  the  yield  and 
ultimate  strengths  of  the  material  in  the  form  of  a safety 
ratio,  SR.  The  safety  ratios  were  comb3jied  by 

I 
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FS- 


5H.f 


(32) 


to  produce  one  safety  factor  for  the  yield  limit  and  one  for 
the  ultimate  limit. 

Results . The  minimum  yield  safety  factor  of  7*^  on 
the  airfoil  occurred  at  mid-sj^an  where  the  maximum  normal 
s tress  was  3298  psi,  and . the  maximum  shear  stress  v/as  2442 
psi.  On  the  strut,  the  minimum  ultimate  safety  factor  v/as 
6.3  at  the  top  of  the  horizontal  portion  v/hich  was  23.973 
in.  from  the  Icaditig  edge  where  the  maximum  normal  stress 
was  9039  psi.  On  the  mounting  boss,  the  mounting  screv/s  had 
the  minimum  ultimate  safety  factor  of  13*3  v/here  the  maximum 
stress  was  4102  psi.  On  the  strut  mounting  plate,  the  mount- 
ing bolts  had  the  minimum  ultimate  safety  factor  of  4.4 
where  the  maximum  normal  stress  was  1032  psi,  and  the  max- 
imum shear  stress  was  7580  psi.  A complete  listing  of  the 
stresses  and  safety  factor's  at  the  critical  points  is  in- 
cluded in  Appendix  C. 

Instrumentation 

Thermocouples,  located  internal  to  the  model,  were  the 
only  instrumentation  selected  for  use  in  the  wind  tunnel. 
Instrumentation  v/hich  was  external  to  the  model  could  not 
be  used  due  to  a wind  tumrel  blockage  problem  tliat  was  ca  ised 
by  the  model  frontal  area.  Additional  internal  instrumenta- 
tion could  not  be  used  since  the  space  inside  tlic  model  was 
very  limited. 
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The  purpose  of  the  thermocouples  was  to  determine  the 
airfoil  surface  temperature  at  selected  points  which  would 
indicate  any  strong  temperature  gradients  that  might  influence 
transition,  establisli  the  cooling  ratio  used  for  that  partic- 
ular test  run,  and  serve  as  a data  base  to  be  used  v/ith 
holography  to  obtain  quantitative  results.  The  thermocouples 
were  needed  primarily  on  the  upper  airfoil  surface  where  the 
boundary  layer  transition  would  be  investigated.  The  follow- 
ing list  of  desired  characteristics  was  devised  to  select 
the  proper  thermocouple  for  this  application: 

1 . Be  very  compact  in  size 

2.  Be  locally  produced  to  meet  specifications 

3.  Useful  range  of  -300  F to  200  F 

4.  Error  range  of  t 25)5  maximum 

5.  Short  time  constant  for  quick  response 

Copper-constantan  thermocouples  were  selected  to  in- 
strument the  model.  This  type  of  thermocouple  was  recommen- 
ded in  Ref  12  for  low  temperature  application  due  to  the 
homogeneity  of  the  metal.  The  thermocouples  were  located  at 
mid-span  of  the  airfoil.  Figure  10  shows  the  location  of  the 
8 thermocouples  along  the  upper  surface . 
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10  Airfoil  Model  Cross  Section 


IV.  Proposed  Test  Procedei-e 


The  wind  tunnel  test  condition 


to  be  used  were  estab 


lished.  Parameters  of  interest  in  tyie  test,  Mach  number 


nolds  .number,  and  cooling  ratio,  were 


;pecified  to  limit  the 


number  of  test  runs  and  to  give  a representative  sample  of 
conditions.  The  Trisonic  v/ind  tunnel  stagnation  temperature 


was  considered  to  be  100  F.  Frees tream  viscosity,  v/hich  is 


a function  of  the  specified  stagnation  temperature  and  Mach 
number,  was  computed  from  Sutherland's  equation  (Eq  10).  Tlie 
only  parameter  left  to  be  computed  was  the  stagnation  pres- 


The  relationship  between  the 


sure 


ubstituting  the 


number,  and  Mach  num.ber  v/; 


nold 


definition  of  Mach  number  and  the  perfect  gas  law  into  the 


definition.  After  introducing  the  isentropic 


Reynolds  numb 


tlie  expression 


ure  and  temperatur 


relationships  for  pre 


ure  wit}i  the  follov/ing 


was  solved  for  the  stagnation  pre 


results 


Equation  33  was  used  to  establish  tlie  proposed  test  condition 
vdiich  appear  in  Table  II . 


Test  Equipment 

Pulse  laser  liolography  and  scliliei'on  were  considered 


Table  II 


Wind  Tunnel  Test  Parameters  lor  R 


n 

M 

(Ib-sec/ft^) 

PQ(psia)  for 
Re-0.923  X 10^ 

P^(psia)  for 
Re-1  .673  X 10^ 

0.3 

0.39279 

9.972 

10.080 

0.5 

0.30358 

6.412 

11 .624 

0.7 

0.37049 

5-058 

9.170 

for  use  during  the  test.  Optical  methods  of  data  collection 
were  necessary  due  to  the  high  blockage  in  the  wind  tunnel. 

The  selection  of  equipment  v/as  based  on  the  information  ob- 
tained from  Ref  10  for  the  schlieren  method  and  Ref  5 for  the 
holographic  method.  Pulse  laser  holography  was  selected  as 
the  primary  method  of  recording  the  boundary  layer  transition. 

The  schlieren  method  would  provide  a recording  of  the 
density  gradients  in  the  flow  field.  Tl:ie  bouiidary  layer  trans- 
ition v/ould  not  appear  as  a point  but  v/ould  appear  as  a 
thickening  of  the  bouiidary  layer  over  a small  region.  Tlie  dis- 
advantages of  tliis  method  were  that  the  results  would  only 
be  qualitative  and  that  the  transition  v/ould  be  difficult  to 
detect  if  t)io  density  gradients  were  not  strong. 

Tlio  pulse  laser  holography  method  liad  several  advan- 
tages. llolngrnpfiy  produces  a hologram  w)iich  cou].d  be  used 
later  to  ^•.'eons  t’-uct  a shadov/grapli , schlieren,  or  intorfero- 

'(•i 


gram.  Qualitative  results  could  be  viewed  in  different  forms 
to  select  the  best  recording  of  transition.  In  additon,  the 
hologram  could  be  magnified  to  help  locate  transition.  An 
interferogram  produced  by  holography  could  be  used  to  pro- 
vide quantitative  results  of  the  density  gradients.  VJith  the 
density  and  tempei-ature  distributions  available,  the  pressure 
distribution  could  be  found. 
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V.  Heat  Transfer  fox’  Tost  Conditions 


There  were  tv.'o  px-oblcms  v/itli  the  desigix  pi'occss  at 
this  point.  First,  the  results  of  Boehman,  Fig  2,  had  not 
been  incox’jxoi'ated  into  the  program,  CONIIEAT.  A critical  Rey- 
nolds number  of  0.5  x 10^,  which  did  not  change  as  a function 
of  the  cooling  ratio,  had  been  used.  Second,  the  pressure 
distribution  used  v/as  for  a Mach  number  of  0.82  since  it 
was  the  only  one  available  for  the  DSMA  523  aii'foil.  Tlie 
model  design  was  completed  by  necessity  with  this  data. 

Test  Data.  Through  a cooperative  test  program  with 
NASA  Ames  Research  Center,  McDonnell  Douglas  Research  Labor- 
atories had  completed  pressure  measux'ements  on  the  DSMA  523 
airfoil  at  a Mach  number  of  0.7  and  Reynolds  number  of  2.0  x 
10^.  Spaid  (Ref  lA)  obtained  special  permission  from  NASA 
and  the  McDonnell  Douglas  Corporation  for  the  data  to  be 
used  in  this  study'-. 

The  new  pressxire  distribution  v/as  ixxcorporated  into 
CONHSAT  along  with  the  vax-'iation  of  the  critical  Reynolds 
number  v^ith  cooling  ratio.  Cai-eful  obsei’vation  of  Fig  2 shov.'S 
that  the  critical  Reynolds  number  may  not  be  reached  fox-  tlie 
given  freestream  Reynolds  nurabei-s  and  cooling  ratios  belov,' 
0.82A.  Therefore,  the  range  of  desired  cooling  i-'atios  v/as 
changed  to  1.000  through  0.82^1. 

Forced  Convection  Fx'edictions . As  a result  of  using 
the  new  test  data  with  CONIIEAT,  predictions  of  'tlie  heat  trans- 
fei’  for  several  cooling  ratios  wei'e  obtained  at  M 0.7»  The 
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variation  of  heat  flux  followed  the  same  pattern  as  the  design 
condition  flux;  hov/ever,  the  magnitude  of  flux  decreased  and 
the  location  of  transition  moved  forv/ard  for  each  increas- 
ing cooling  ratio.  As  can  he  seen  in  Fig  11,  the  transition 
to  turbulent  flow  did  not  occur  at  all  for  Re=0.923  x 1 0^ 
and  Ty/i'co  -0.824.  The  remaining  heat  flux  predictions  can  be 
found  in  Appendix  B. 

Due  to  the  differences  in  design  and  test  data,  a com- 
parison of  the  total  heat  transfer  rate  for  these  two  condi- 
tions was  made  to  ensure  the  capacity  of  the  cooling  system 
was  adequate  for  the  test.  This  comparison  (Fig  12)  showed 
that  the  required  cooling  for  both  freostream  Reynolds  num- 
bers was  less  than  the  design  values  at  the  same  cooling 
ratio . 
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flZRFOlL  HEAT  FLUX  I BTU/SEC-SQF ) 


X/C  ILOCfUlON/CHORG  J 


Fig.  11  Airfoil  Heat  Flux  at  M=0.7,  Re=0.923  x 10^,  and 
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COOLING  RATIO  (T ) 


1 


The  follov/ing  items  summarize  the  capability  of  the 


airfoil  test  model: 

1.  The  angle  of  attack  v/hich  is  not  adjustable 
is  ^ = 0°. 

2.  The  m.aximum  allowable  aerodynamic  loads  occur  at 
M=0.7. 

3.  The  maximum  cooling  capacity  is  17.03  B/sec. 

4.  The  maximum  moss  flow  r-ate  of  coolant  is 

m=0.198  lb  /sec. 
m' 

5.  The  maximum  temperature  gradient  in  the  model 

is  2.23  F. 

6.  The  maximum  cooling  experienced  by  the  wind  tunnel 
window  is  0.004  B/sec. 

7.  The  lowest  safety  factor  v/hich  occurs  at  tlie  strut 
mounting  bolts  is  4.4. 

The  first  four  recommendations  are  concerned  with  the 
testing  of  the  airfoil  model.  The  remaining  recommendations 
propose  further  investigations.  These  recommendations  are: 

1.  A bench  check  of  the  model  cooling  system  and  ther- 
mocouples should  be  made  prior  to  the  test. 

2.  Ensure  the  wind  tunnel  is  dry  as  possible  to  pre- 
vent frost  formation. 

3.  Compare  the  total  heat  flow  for  the  experiment  with 


predicted  values  by  measuring  the  mass  flow  rate  of 
coolant. 


4.  Compare  the  experimental  results  with  the  theoret- 
ical results  from  Ref  2. 

5.  Further  experimental  investigation  of  the  cooling 
effects  on  transition  should  be  conducted  v/ith  a 
flat  plate  to  eliminate  pressure  gradients. 

.6.  Further  experimental  investigation  of  the  cooling 


effects  on  transition  should  be  conducted  in  a 
three-dimensional  flov/. 
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Appendix  A 

Computer  Tro^ram  for  Boundary  Layer 
Convection  Study 


The  computer  program  called  CONHEAT  implemented  the 
theory;  and  program  logic  that  were  discussed  in  Section  II 
to  perform  the  convective  heat  transfer  analysis.  To  use 
this  program  effectively,  the  format  for  input  and  output 
and  definitions  of  parameters  must  be  Inov/n.  A listing  of 
the  complete  program  has  been  included  in  this  appendix. 

Input  Format  and  Parameters 

Two  methods  wore  used  to  input  data.  One  method  was 
the  data  statement,  and  the  other  was  standard  data  cards. 

The  data  statement  is  part  of  the  main  progr-am  and  was  used 
for  constants  and  variables  which  wore  seldom  changed.  Data 
cards  v/ere  used  for  information  v/hich  was  generally  changed 
or  for  variables  which  had  to  be  redefined  during  a single 
execution  of  the  program. 

Five  data  statements  were  used.  The  data  statement 
assigns  values  to  variables  and  arrays.  The  v/ord  "DATA"  may 
be  followed  by  any  number  of  variables  with  the  value  of  each 
enclosed  with  slash  marks.  The  first  data  statement  defined 
constants  and  geometry  of  the  model  as  follows: 

GO  - Newton's  constant,  gQ=32.17^  Ib^^^-f t/lb^-sec^ 

R - Gas  constant  for  air,  R=17l6  ft^/sec^-R 
CJ  - Mechanical  to  thermal  energy  conversion  factor, 
J-778.16  ft-lb^B 
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C3H-  opecific  heat  at  constant  preacure,  0^=0. 24 
PR  - Prandtl  number  for  air,  Pr=0.7 
C - Airfoil  chord  in  feet,  c=0. 75  ft 
B - Airfoil  span  in  feet,  b=1.83  ft 
RA  - Leading  edge  radius  in  feet,  r=0.025  ft 
TIN-  Freestream  temperature,  -559*67  R 
The  second  data  statement  defined  the  distance  between  pres- 
sure coefficient  locations  on  t}ie  upper  surface  in  the  sur- 
face coordinate  system.  XU  means  the  x-distance  on  the  upper- 
surface.  All  21  points  must  be  included.  The  values  in  this 
array  are  separated  by  commas.  The  follov/ing  is  an  example: 

DATAXU/0.,  .048,  .086,  ...,  .7?0/ 

The  third  data  statement  has  the  same  form  as  the  second 
statement.  It  defines  the  x-distance  on  the  lower  surface, 

XL.  The  fourth  data  statement  lists  the  pressure  coefficients 
for  the  upper  surface,  CPU.  The  form  is  the  same  as  the  sec- 
ond statement.  The  fifth  data  statement  lists  the  pressure 
coefficients  on  the  lov/er  surface,  CPL . 

The  three  read  statements  in  the  program  are  unformated 
to  simplify  use.  This  requires  the  data  card  to  list  the 
decimal  values  of  the  variables  separated  by  commas  and  in 
the  order  that  appears  on  the  read  statement.  The  first  data 
card  must  contain  the  following: 

UIN  - Freestream  velocity  in  ft/sec,  u^^ 

MUIN-  Freestream  viscosity  in  lb^,-sec/f  t^ 

PIN  - Freestream  pressure  in  Ib^/ft',  P^^, 
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An  example  of  this  card  is  the  following: 


813.052,  .00000037049,  728.048 
Tiie  numhers  start  in  column  1 of  the  card.  Thex-e  must  be 
at  least  one  card  of  this  type.  The  second  card  indicates 
the  number  of  cooling  ratios  to  be  read,  NTR.  The  value  on 
the  second  card  must  bo  a positive  integer.  The  third  data 
card  lists  the  following: 

TH  - Cooling  ratio,  T^T^^ 

REG  - Critical  Reynolds  number  corresponding  to  the 
value  of  TR  of  Fig  2 

The  total  number  of  the  third  tj^e  of  cards  must  equal  the 
integer  on  the  second  card.  If  more  than  one  set  of  fx-ee- 
stream  conditions  is  specified,  then  all  three  types  of 
cai’ds  must  bo  repeated. 

0utT'>ut  Format 

Two  types  of  output  ai-'e  obtained  from  the  program.  One 
type  is  a plot  of  the  local  heat  flux,  q^,  as  a function  of 
location  al.ong  the  chox'dline  for  the  upper  and  lov/er  sur- 
faces. ExEimples  of  the  plot  appear  in  Appendix  B.  The  second 
type  of  output  is  a Uiblc . The  table  includes  freestream 
conditions,  local  values  of  propex'ties,  and  the  total  i-ate 
of  heat  transfer.  The  local  values  of  propex’ties  at  each 
point  Include  location,  velocity,  Reynolds  number,  Stanton 
number,  convective  conductance,  and  heat  flux.  An  example  of 
this  table  appears  in  Appendix  B,  Table  III. 
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Appendix  B 

Convection  Prof^ram  Output 

The  computer  program,  CONilEAT , which  v/as  discussed  in 
Section  II  and  Appendix  A v/as  run  with  the  proposed  test 
data  discussed  in  Section  V.  These  results  predict  the  heat 
transfer  and  flow  properties  for  the  proposed  test  conditions 
at  a Mach  number  of  0.7.  An  example  of  the  tabular  results 
is  contained  in  Table  III.  Plots  of  heat  flux  for  all  of 
the  proposed  cooling  ratios  are  contained  in  Fig  13  through 
Fig  19. 
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Table  III 


Airfoil  Local  Properties  at  the  Test  Condition 
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O UPPER  SURFACE 
A LOWER  SURFACE 


Fig.l4  Airfoil  Heat  Flux  at  M=-0.?,  Re -0.923  x 1 and 
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Fig.  15  Airfoil  Heat  Flux  at  IVI-O.7,  Re  = 0.923  x 10  , and 
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1 1 1 1 7 I 1 r 1 1 

00  O.IO  O.ZO  0.30  0.40  O.CO  O.CO  0.70  0.60  O.SO  1.00 


I 


O UPPER  SURFnCE 
A LOWER  SURFRCE 


Fig.  16  Airfoil  Heat  Flux  at  M=0.7,  Re=l .673  x 10^,  and 
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ZRrOIL  HEAT  FLUX  (6TU/SEC 
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flZRfOZL  Hfar  FLUX  {aTU/SEC-SQD 


X/C  ILOCfiTION/CHORD) 


Fig.  19  Aii-'foil  Heat  Flux  at  M=--0.7,  Re .673  x 10^,  and 
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Stress  Analysis  Cajculations 

Airfoil 

IJoments  of  Inert.]  a and  Centroids . With  the  aid  of  a 
full  scale  drawing  of  the  airfoil  section  such  as  Fig  10, 
the  area,  centroid,  and  moments  of  inertia  about  the  leading 
edge  v/ere  computed  using  Eq  21-26.  Then,  the  parallel  axis 
theorem  was  applied  to  get  the  moments  of  inertia  for  axes 
through  the  centroid  and  for  the  x-axis  parallel  to  the  chord. 
The  results  v/ere: 

A = 3.^30  in.^ 

= 3.9^^7  in. 

C = 0.042  in. 

J 

I = 0.118?  in.^ 

I =17.0948  in.^ 

I = 0.3940  in.^ 

material  Selection.  The  inate’'ialc  used  for  tlie  con- 
struction of  the  airfoil  and  end  sections  are  listed  in  Table 
VII.  A complete  listing  of  structui’al  properties  needed  to 
perform  the  entire  stress  analj^sis  is  included  in  Table  IV 
ty  'type  of  material . 

Equivalent  Force  System.  The  computed  uniform  lift 
force  at  the  design  condition  was  12.458  Ib^./in . of  span  for 
a total  of  299  Ib^.  normal  to  the  chord.  The  uniformly  distrib- 
uted drag  force  was  0.537  Ib^/in . of  span  for  a total  force 
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Table  IV 


Structural  Limits  of  Materials 


Type  Material 

•dr 

Tension 

1 

C ompre 

* 

ssion 

* I 

Shear 

Yld  Ult 

Yld 

Ult 

Yld 

Ult  I 

1 

202^'-T4  Alum 
(Ref  -9:85) 

4o  64 

38 

@ 

23 

4o 

AISI-SAE-4130 
SteeKRcf  14:452) 

117  136 

@ 

136 

© 

102 

Oak,  ViJhite 
(Ref  9:146) 

@ 0.8 

@ 

@ 

@ 

2 

Yld  - Yield 

All  stresses  in 

1000  psi 

Ult  - Ultimate  @ Not  available 

of  12.9  Ib^.  The  moment  of  672  In.-lb^  was  distributed  uni- 
formly at  28.00  in.lb^/in.  of  span.  These  are  tiie  tlieoretical 
loads  computed  in  Section  III.  This  system  of  forces  was 
concuri'ent  at  the  quarter  cliord  point. 

Critical  Points.  Due  to  model  symmetry,  only  half  of 
the  airfoil  had  to  be  analysed.  Thei'c*  w'ore  5 critical  points 
selected; 

1 . Mid -span 

2.  4.5  in.  from  mid-span 

3.  Center  of  airfoil  boss 

4.  11.0  in.  from  mid-span 

5 . End  of  span 

Special  Considerations.  The  elliptical  shell  which 
replaced  the  airfoil  section  for  the  analysis  b.ad  a wall 
thickness  of  0.229  in.  for  a 9 in.  major  axis  and  0.99  in. 
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minor  axia. 


The  ellipse  was  defined  by  solving  the  follov/ing 


equation  for  thickness,  t: 

Area  of  Ellipse  — [ ah’^C y'X.'Jj  ~ 3/i 3 in/ 

Equation  27  was  used  to  compute  the  sheai'  stress  at  critical 
points  1 through  ^4-.  Equation  28  was  used  at  point  number  5 
for  the  shear  stress  of  the  solid  elliptical  section. 

Stresses  and  Safety  Factors.  The  airfoil  was  assumed 
to  behave  as  a simple  beam  with  a uniform  load  and  symmetri- 
cally overhanging  ends.  The  normal  and  shear  stresses  were 
computed  for  4 locations  at  each  critical  points  the  leading 
edge,  trailing  edge,  top  of  the  section  at  maximum  thickness, 
and  bottom  of  tlic  section  at  maximum  thickness.  The  stresses 
were  combined  using  Eq  30  and  3i ■ The  safety  factors  were 
computed  v;ith  Eq  32.  These  results  are  in  Table  V. 

Mounting  Struts 

Moments  of  Inertia  and  Centroids.  The  vertical  portion 

of  the  strut  was  assumed  to  resist  bending  v/ith  the  exhaust 

pipe  0]ily.  The  moment  of  inertia  of  tlie  3/^  in.,  nominal 

sized  pipe  was  0.0448  in.^  The  moment  of  inertia  niid  centroid 

of  the  horizontal  portion  was  computed  using  Eq  21-26  with 

the  drawing  of  the  cross-sectioi  of  tlie  strut  in  Fig  7»  Tlie 

result  was  an  area  of  0.8^402  in.  , the  centroid  vhiich  was 

0.500  in.  from  the  center  of  tlie  exhaust  piT:)e,  and  the  mo- 

4 

ment  of  inertia  wliich  was  0.272  in.  about  tlie  axis  perpen- 
dicular to  the  axis  of  symmetry  of  the  section. 

Material  Selection.  The  materials  used  for  the  con- 
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Table  V 


Maximum 

Stresses 

and  Safety 

Factors 

— 

for  the 

Airfoil 

Critical  Point 

x--r^ 

max 

^max 

F3 

FS  i 

1 

1 

LOCATION 

lO^psi 

p 

lO'psi 

YIELD 

ULTir.iATE 

Top 

32.973 

^24.4l8 

7.4 

1 

12.5  ! 

1 

Bottom 

14.873 

25.101 

8.7 

@ 1 

Front 

, 0.461 

0.483 

4-15.2 

@ 

Back 

0.512 

0.483 

408.3 

1 

707.1 

Top 

19.140 

1 2 . 991 

13.5 

22.7  ' 

2 

Bottom 

6.279 

13.571 

16.3 

@ 

Front 

0.389 

0.405 

500.0 

Back 

0.426 

0.405 

488.0 

845.2  1 

Top 

0.086 

1 .155 

200.0 

1 

3 

Bottom 

2.609 

1 .431 

114.4 

189.7  1 

1 

Front 

0.439 

0.436 

46l  .3 

816.5 

Back 

0.432 

0.436 

466 . 3 

1 

0 

1 

Tox' 

4.967 

5. 086 

39.4 

1 

@ 1 

4 

Bottom 

5.228 

5. 08? 

38.1 

66.2  i 

1 

Front 

0.144 

0.144 

1581 .3 

2786.3  1 

Back 

0.143 

0.144 

1581  .3 

5 

Top 

0. 

0.485 

@ 

4l  .2 

@ Not  computed 


I 


I 


i 

i 
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struction  of  tlie  strut  are  listed  in  Table  Vli.  The  structui- 
al  limits  of  these  materials  are  in  Table  IV.  The  epoxy  eas- 
ing material  was  not  used  to  support  loads  in  the  analysis. 

Equivalent  Force  System.  Duo  to  the  model  syinmotry, 
only  }ialf  of  the  loads  used  to  analyze  the  airfoil  acted  on 
each  vertical  portion  of  t)ie  strut.  Ttie  strut  drag  of  7*27 
Ib^.  which  acted  at  the  middle  of  the  strut  was  added  to  the 
system.  To  obtain  a concuirent  foi-ce  system,  the  strut  drag 
v/as  moved  parallel  to  the  location  of  the  airfoil  ctiordline, 
and  a couple  equal  to  the  drag  force  times  ttie  translated 
distance  v/as  also  added.  The  resulting  equivilent  system 
was  the  follov/ing; 

L ^ Ibj. 

D - 13.7  Ib^ 

M 291  .5  in.-lb^. 

This  equivalent  system  for  the  vpi'tical  portion  was  traiiS- 
lated  in  a similar  manner  to  obtain  the  foi’cus  an>’  inorm nt 
acting  on  tlie  horizontal  i)oi’ticjii.  U]ion  transla tita  to  the 
base  of  the  vertical  portion,  tJie  rosuitinf,  ((paivulent  sys- 
tem for  the  horizontal  portion  v/as: 

L - 149.5  ibj. 

D = 13.7  Ib^ 

M = 442.4  in.-lb^ 

Critical  Points.  There  were  only  two  critical  points 
to  analyze  on  the  strut.  These  were  the  points  of  maximum 
bending  moment.  One  point  was  11.0  in.  down  from  the  airfoil 
chordline  on  the  vertical  portion  of  tlie  strut.  Tlio  otlror 
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point  v/aa  P-')>97^  in.  frani  the  C(./»lor  of  tho  verti'.-al  i.xJjaual 
pipe  . 

Special  Jon"  i't»  i - ^ I Lon;; . The.  conputci  norruil  on 

the  ver  tical  portion  of  the  .strut  wajj  eorr-  oteJ  for’  cm  vnt 
affectn.  The  nornai  ntt*j;}  viCi.;  n ul  tirl  i*ri  fy  a con '.•e  t i i r 
factor  on  tho  iriside  an  l (uit  1 tlic  f wnioi.  An  t-  l .7i. 

and  0.71.  r»  tlv»  ly  . 
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mum  normal  sUi  su  was 
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1 r mar.  ■ ■ . 

11 

• . ! a; 

was  4275  p3 i . anJ  the 

' 1 i i 

- i* 

...  ‘‘a  • . 

tr'-.n  1 

) . t ( 

the  oi'it.siile  of  the  r 1' 

1 nti  . 

, t r,  I)  -., 

l.or 

f 1 a t ■ « ; I 

was  3916  pa  i , th-’  max 

i 'uj’ 

.rhea-  : 

t ; a;  w'-..  ' 

■ ’ ?'-J 

1,  .M  l the 

ultiiiiato  naff'ty  f.ictoi  w.:  28.  ^■9.  si  tti<  t ;(  f tfr  I..  i ' ri- 
tal  strut,  tho  tnaximum  normal  :.ir.  wau  i ,0,6  T..i,  Ih- 
maximum  shear  stress  was  9030  js  ’ , onJ  th*'  ultiwat  - i;.»''.-*y 
factor  was  6.'J.  On  t)a'  I'ottom  of  t)i"  hor  iv.on ta]  strut,  the 
maximum  normal  stress  v.as  i-si,  the  ^.a.^iinum  shear  rttoss 

was  1580  pai,  and  the  safi;  ty  factor  v/as 

Mounting;  Boos 

— - - - --i 

Moments  of  Inertia  and  Centroids.  Duo  to  tlio  nature  of 
the  assumed  equivalent  force  system,  no  moments  of  inertia 
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or  ccnti-oiclu  wt  i-',*  calcula  U d . 


Mut<!riul  ;:<»Ii‘ctiori . The  matfrial  uoed  to  coni' true  t 
the  n.ouriting  bon.;  wan  20?h  T4  aluminuir.  Tlii:;  material':: 
stru'^tural  limit,.:  are  lisite-l  in  Tut  li  IV.  The  6-'f0  nKiclilne 
acr«*wj  'ii;od  ha.l  jw;  ultinu.te  utrnoj  limit  oT  yi,400  pni.  The 
ultir..;l"  Btieai  utioss  foi"  the  fill"!  v/olci  wa::  11,300  j ::  i . 

lv'<l  e:;  t ’I' J>.  ' V m . Tha  •*amo  li-u  1j  l.hat  ectud 
on  tha  viti  I tortion  it  llie  stiut  aioo  U't.l  on  Mu  inoun- 


ting  bv< 

tiH  . H'lwev 

1 , ttic  d:  Fg  was 

/ii* 

elected  si) 

CO  it  w.n. 
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. Tte  ’ •-  w^i  !■ 

three  riti 

eai  j'oints 

nn 

1 ^ ‘ ri  ■ . 

t.iWu  b V 

. rheo"  t'l'intu 

wtf : 

tt.i 

which 

i' t- 

1 ir 

fT  ^ i ■ 

* !ie  n i r f . n PnJ 

boss  plnlu, 

•i"  rimigi 

ef 

tfK>  1.  a:;,  and  tf  a 1 -nF,  jlBte  flllot  V.  1 i . 

oltewi^^n  fin  ! .laT^ty  ■'•  Ic.-.t  . I’horo  were  11  fjcrewc 

th"  butn;  to  th'  ail  foil  5 11  ncr-.  far?  lur  the 

bouB  to  ttf>  tom  plate.  The  mnAimur  Btien.a  i.;i  the  fii’rowj  ivah 
^♦102  p::  , and  tin*  ullimate  nafety  factor  wan  13*3«  'fl*'  maxi- 
mum ohrux-  atrear;  oji  ttie  flan^^i*  of  the  bona  wliicli  waa  correct- 
ed for  ntrenn  concentration  was  890  psi.  The  yield  safety  fac- 
tor was  44.9.  The  maximum  shear  otres.s  on  the  fillet  weld 
was  1032  psi,  and  the  ultimate  safety  factor  wa.s  43.8.  These 
strGS.ses  and  safety  factors  wore  computed  wit)i  Eq  .?9  ttiroupih 
Eq  32. 
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l>*ountinK  Fla  be 

Mon.nnts  of  Inertia  and  Ci-introidr; . Due  to  the  nature  of 
the  equivalent  force  uyaten,  no  moments  of  inertia  or  cen- 
troids were  calculated. 

Materiel  Selection.  The  material  used  to  construct 
the  mounting  plat(j  wan  4l30  low  alloy  steel.  This  material's 
sLructural  limits  are  listed  in  Table  IV.  The  fillet  weld  had 
an  ultimate  shear  limit  of  43,200  jisi.  Tiie  mountin/r;  bolts, 

1/4  in.  hexiiead,  had  an  ultimate  normal  stress  of  4^1 , 800  psi 
and  an  ultimate  shear  stress  of  'J3,600  j'si. 

Equivalent  Force  System.  The  equivalent  force  system 
was  the  same  as  the  loads  that  acted  on  the  liori^.ontal  })Or- 
tion  of  the  strut. 

Grit' cal  Points.  There  were  two  critical  jioints.  These 
points  weir  tiie  l/8  in.  fillet  wold  between  the  plate  an;I 
thv  strut  and  the  mounting  bolts. 

Strosnc!-.  and  Safety  Factors.  The  s trusses  and  safety 
factors  v/ei-o  computed  v/ith  Eq  29  through  Eq  32.  The  maximum 
shear  otrers  on  the  fillet  v/cld  was  VOSO  psI,  and  the  ulti- 
mate safety  factor  v.’as  6.0.  The  maximum  normal  stress  on  t]\o 
mounting  bolts  was  1 032  psi,  and  the  maximum  slicar  stress  v/as 
7580  psi.  The  ultimate  safety  factor  on  the  bolts  \/;is  4.4. 
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Appendix  D 


Details  of  Model  Construction 

As  a culmination  of  the  design  and  analysis  process, 
details  necessary  for  the  construction  of  the  model  have 
been  specified.  P’igure  5.  7,  8,  9i  10|  20,  and  21  depict  all 
of  the  major  parts  that  comprise  the  airfoil  and  strut  mount 
system.  Only  Fig  20  and  2l  appear  in  this  appendix.  Figure 
20  contains  a perspective  of  the  entire  assembly.  Notes  which 
explain  special  cojisideratioriS  during  construction  or  clarify 
details  }iavo  been  included  in  tlie  drawings.  The  airfoil  sec- 
tion shape  used  is  the  McDonnel.!  Douglas  Corporation,  DSMA 
523,  su]^ercritical  airfoil.  Table  VT  provides  geometric  coor- 
dinates (Ref  8:738)  necessary  to  contour  tl^e  surface  of  the 
airfoil  for  a 9 inch  chord. 

Sevei-al  itcr'.s  require  special  attention  during  construc- 
tion. Teflon  gaskets  must  be  used  on  all  contact  surfaces 
between  parts  to  pi-cvent  leakage  of  liquid  and  gaseous  nitro- 
gen. After  tlic  thermocouples  have  been  installed  using  the 
1 0-32  screw.s  wit}i  the  center's  drilled  as  seen  iji  Detail  A of 
Fig  10,  the  centers  must  be  back  filled  with  epoxy  to  secure 
tlie  tliermocouple  wii'cs  and  insulate  the-  thermocouple  from 
direct  contact  wit]i  the  liquid  nitrogen. 

Lastly,  aluminuraized  milar  and  fiber  glass  multilayered 
insulation  must  be  used  to  v;rap  all  thermocouple  wires  and 
the  liquid  nitrogen  supply  line  that  iedidG  to  the  mounting 
strut . 
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Table  VI 


Airfoil  Section  Coordinates 


X 

t 

X 

t 

.0045 

.0456 

-.0459 

.0915  ' * 

4 .i4oo 

.4963 

-.4693 

.9656 

.0090 

.0639 

- .0642 

.1281  ■ ! 

4.3200 

.4948 

-.4602 

.9550 

.0225 

.0996 

-.0997 

.1993  ! i 

1 1 

4.5000 

.4925 

- .4492 

.9417 

.0450 

.1379 

-.1379 

.2758 

1 . 

i 4.6800 

.4895 

-.4363 

.9258 

i .0675 

.1658 

-.1658 

.3316 

; 4.8600 

.4858 

-.4210 

.9068 

.0900 

.1864 

-.I860 

.3724 

j 5.0400 

.4814 

-.4015 

.8829 

.1125 

.2039 

- .2029 

.4068 

I 

5.2200 

.4763 

-.3781 

.8544 

.1350 

.21 84 

-.2172 

.4356 

5 .4000 

.470'-,' 

-.3500 

.8204 

.1800 

.2423 

-.2407 

.4830 

5.5800 

.4639 

-.3166 

.7805 

.2700 

.2766 

-.2760 

.5526 

5.7600 

.4565 

-.2785 

.7350 

.3600 

.3011 

-.3025 

.6036 

5.9400 

.4484 

-.2375 

.6859 

.5400 

.3367 

-.3428 

.6795 

6 .1200 

.4395 

-.1939 

.6334 

.7200 

.3633 

-.3757 

.7390 

6.3000 

.4295 

-.1526 

.5821 

1 .9000 

.3869 

-.4009 

.7878  , 

6.4800 

.41 84 

-.1142 

.5326  i 

1 1 .0800 

.4068 

-.4212 

.8280 

6 .6600 

.4062 

-.0788 

.4850 

! 1 .2600 

.4232 

-.4375 

.8607 

6.8400 

.3926 

- . 046  8 

.4394 

j 1 .4400 

.4369 

-.4510 

.8879 

7.0200 

.3775 

- . 01 84 

.3959 

1 1 .6200 

.4485 

-.4621 

.9106 

7.2000 

.3604 

.0062 

.3542 

1 1 .8000 

.4581 

-.4713 

.9294 

7.3800 

. 341 2 

.0267 

.3145 

: 1 .9000 

.4662 

-.4789 

.9^51  i 

7.5600 

.3195 

.0428 

.2767 

1 2.1600 

.4731 

-.4850 

• 9581 

7.7400 

.2950 

.0542 

.2408 

2.C400 

.4?88 

- .4808 

. 9686 

7.9200 

.2670 

.0602 

.2068 

2.5200 

.4836 

-.4933 

.9769 

8.1000 

.235^ 

. 0595 

.1759 

2.7000 

.4874 

-.4955 

.9829 

8.2800 

.1997 

.0507 

.1490 

2.8800 

.4906 

-.4965 

.9871 

8.4600 

.159^ 

.0321 

.1273 

3.0600 

.4931 

-.4962 

.9893 

8.6400 

.1138 

.0031 

.1107 

3.2400 

.4950 

-.4948 

.9898 

8.8200 

.0616 

-.0379 

.0995 

3.4200 

.4963 

-.4923 

.9886 

9.0000 

.0028 

-.0910 

.0938 

3.6000 

.4972 

-.4885 

.9857 

3.7800 

.4974 

-.4835 

.9809 

L 



All  measvirements  in  inches  (accurate  to  0.0001  in.) 
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Mater ialLi  to  be  used  for  construction  of  the  model  were 
determined  from  consideration  of  the  heat  transfer  and  stress 
analysis.  The  thermal  and  structural  properties  of  these 
materials  v/ere  required  to  exceed  the  limits  of  the  design 
condition.  Table  VII  lists  materials  to  bo  used  for  the  com.- 
ponents  of  the  model. 


Table  VII 


Specification  of  Construction  Materials 


Component 

Airfoil 

Center  Section 
End  Plates 

Insulating  End  Sections 
Airfoil  Bosses 
Cooling  Nozzles 
Strut  Boss  Plates 
Mounting  Struts 


Material 

2024-T4  Aluminum 

2024-T4  Aluminum 

VJhite  Oak  (Hardwood) 

2024-T4  Aluminum 

Type  30^+  Stainless  Steel 

AISI-SAE4130  Low  Alloy 
Steel 


Exhaust  Pipes 
Thermocouple  and  LNg  Pipe 
Reinforcing  Bars 
Contour  Shape 
Strut  Mounting  Plates 
Thermocouple  and  LN2  Supply  Pipe 


AISI~SAE4130  Low  Alloy 
Steel 

AISI-SAE'^H30  Low  Alloy 
Steel 

AISI-SAE^130  Low  Alloy 
Steel 

C30I  Epoxy  (Aluminum 
Filled) 

AISI-SAE  Low  Alloy 
Steel 

Type  30^J-  Stainless  Steel 


The  requirem.ents  for  fastening  methods  wore  based  on 


/] 
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analysis . 


the  stress 


Tv/o  primary  me  tliods  which  will  be  em- 


ployed are  fillister  head,  fine  thread,  machine  screv/s  and 
fillet  welds.  Table  VIII  lists  the  fastening  method  with  sizes 
for  all  of  the  components  of  the  model. 


Table  VIII 

Specification  of  Fasteners 

Component  Combination  T^goe  Fastener  Size 

Airfoil  Center  Section  Ka.lves  Fillister  Screv/  10-32 

Airfoil  End  Plates  and  Insulating 

End  Sections  to  Center  Section  Fillister  Screw  10-32 

Airfoil  Boss  to  Airfoil  Center 

Section  Fillister  Screw  6-40 

Cooling  Nozzles  to  Airfoil  Boss  Fillister  Screw  6-40 

Airfoil  Boss  to  Strut  Boss  Plate  Fillister  Screv/  6-40 

Strut  Boss  Plato  to  Strut  Pipe  Fillet  VJeld  1/8  in. 

Reinfox’cing  Ears  to  Exhaust  Pipe  Fillet  Weld  3/l6  in. 

Reinforcing  Ears  to  LN2  Pipe  Fillet  V/eld  3/32  in. 

Struts  to  Mounting  Plate  Fillet  VJeld  1/8  in. 

Mounting  Plate  to  Wind  Tunnel  Hexhead  Screw  1/4  in. 

28  Threads 
per  inch 
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^A  wind  tunnel  test  model  of  a supei’critioal  airioil  was 
designed  to  investigate  the  wall  cooling  effect  on  subsonic 
boundary  layer  stability.  A DSMA  523  airfoil  section  was  em- 
ployed. The  model  w.as  desired  to  have  surface  temperature 
instrumentation  and  a liquid  nitrogen  cooling  system. 

Hoat  transfer,  aerodynamic  loads  and  stresses,  and 
instrumentation  were  analysed  for  the  pro])osed  test  condi- 
tions. A computer  program  was  developed  to  analyse  the  forced, 
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-T^convcctivc  boat  transfer  over  a two-dimensional  body  v/ith  a 
constant  wall  temperature.  The  program  utilized  an  integral 
motliod  to  compute  local  Stantoti  numbers.  Local  heat  flux  and 
total  heat  flow  v/ere  predicted  for  a Mach  number  of  0.7i  Hey- 
nolds  numbers  of  0.9^3  x {109  and  1 .673  x and  cooling 

ratios  from  \ .000  to  0.8247  The  stress  analysis  consisted  of 
applying  beam  bending  theory,  along  with  some  simplifying 
assumptions,  to  the  model.  Construction  drawings  and  speci- 
fied test  conditions  for  Mach  numbers  of  O.3*  0.5,  and  O.7 
are  included.  The  proposed  test  are  to  be  conducted  in  the 
subsonic  test  section  of  the  Trisonic  Test  P'acility  at  VJrigiit 
Patterson  AFB,  Ohio. 
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